Three-body loss in lithium from functional renormalization 
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We use functional integral methods for an estimate of the three-body loss in a three-component 
6 Li ultracold atom gas. We advocate a simple picture where the loss proceeds by the formation of 
a three-atom bound state, the trion. In turn, the effective amplitude for the trion formation from 
three atoms is estimated from a simple effective boson exchange process. The energy gap of the 
trion and other key quantities for the loss coefficient are computed in a functional renormalization 
group framework. 



Ultracold fermion gases with three components show 
new features as compared to the well studies systems with 
two components. For degenerate fermions with SU(3) 
symmetry one finds in the unitarity limit of infinite scat- 
tering length the interesting tower of Efimov states []]]. 
These are series of three-atom bound states with a geo- 
metrically decreasing gap parameter, reflecting the vio- 
lation of scale symmetry by a limit cycle scaling behav- 
ior of the renormalization flow @, Q ■ In the vicinity of 
the Feshbach resonance the lowest three-atom state, the 
trion, is below the open channel energy level. At low tem- 
perature one can infer a phase structure different from 
the BCS-BEC crossover for the two-component system, 
namely an intermediate trion phase without superfluid- 
ity separating the superfiuid BCS and BEC phases 
Interesting quantum phase transitions may describe the 
phase transition between phases at vanishing tempera- 
ture. 

The trion bound state is also expected to persist if the 
SU(3) symmetry is violated by a different location and 
strength for the Feshbach resonances between different 
pairs of atomic components. Recent measurements of 
the three-body loss coefficient in a three-component sys- 
tem of 6 Li 0, H| may find an interpretation in this way 
0, 0|- We investigate here a simple setting, where the 
loss arises from the formation of an intermediate trion 
bound state, which subsequently decays into unspecified 
degrees of freedom - possibly the "molecule type" dimers 
associated to the nearby Feshbach resonances. In turn, 
the trion formation from three atoms proceeds by the ex- 
change of an effective bosonic field, as shown in Fig. [1] 
We estimate the loss coefficient K$ as being proportional 
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Here m x and T x are the trion gap parameter and decay 
width, while describes a type of gap parameter for 
the effective boson, such that its propagator can be ap- 
proximated by rn^~ . The Yukawa couplings hi couple 
the fermionic atoms to the effective boson, and the trion 
coupling <7i accounts for the coupling between trion, atom 
and effective boson. We sum over the "flavor" indices 
i = 1,2,3. We will estimate m x , nrL, hi and gi from the 
non-perturbative renormalization flow which arises from 
a simple truncation of the exact flow equation for the 
average action or flowing action 0], for reviews see 

Recently we used the method of functional renormal- 
ization to describe a SU(3) invariant system of three 
fermion species close to a common Feshbach resonance 
[H, floj | . In this context we explored the manifestation 
of the Efimov effect and formulated some predictions on 
the quantum phase diagram in such systems. In contrast 
to this theoretical model, the system consisting of three- 
component 6 Li atoms, which is of current experimental 
interest [3, @|, does not possess this SU(3) symmetry. 
The main difference is that the resonances do not occur 
at the same magnetic field, and thus, for a given mag- 
netic field B, the scattering lengths of different pairs of 
atoms, (1, 2), (2, 3), and (3, 1) differ from each other. 

In this letter we generalize the model presented in Q 
to cope with this more general situation. Our truncation 
of the (euclidean) average action then reads 
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FIG. 1: Three-body loss process involving the trion. 
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where we choose natural units h — 2M — 1, with the 
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atom mass M. We sum over the indices i, j, k wherever 
they appear. Here ipi denotes the fermionic atoms, 4n a 
bosonic auxiliary field which mediates the four-fermion 
interaction and \ is a fermionic field representing the 
bound state of three atoms. Formally, this trion field 
is introduced as the field mediating the interaction be- 
tween atoms ip and bosons <p>. We show this schemat- 
ically in Fig. [2] In the limit m 2 — > oo, to 2 ^ — > oo, 

h-i/ m % ~~ * l-M> 9i/ m x ~ * I'^l tne ac tion describes 
pointlike two-body interactions with strength Aj, as well 
as a three-body interaction with strength A^. We will 
concentrate on a microscopic interaction of this pointlike 
type. We consider the "vacuum limit" where tempera- 
ture and atom density go to zero. Then the chemical 
potential /i in Eq. |2]| satisfies /i < 0. A negative chemi- 
cal potential \i has the meaning of an energy gap for the 
fermions when some other particle (boson or trion) has 
a lower energy. The dominant difference to the SU(3) 
symmetric model arises from the different propagators of 
the bosonic fields 4>i=ip2ip3i </ , 2=V'3'0ij and 4>z=ip\ip2- In 
addition, we allow in general for different Yukawa cou- 
plings hi corresponding to different widths of the three 
resonances. Also the Yukawa-like coupling gi that cou- 
ples the different combinations of fermions ipi and bosons 
4>i to the trion field x=4'i' t p2'4'3 is permitted to vary with 
the species involved. Although the SU(3) symmetry is 
explicitly broken, the system exhibits three global U(l) 
symmetries corresponding to the three conserved num- 
bers of species of atoms. 





FIG. 2: Interaction between atoms ip and effective bosons <p 
as mediated by the trion field \- 

The rcnormalization flow of the various couplings from 
the microscopic (UV), k = A, to the physical, macro- 
scopic (IR) scale, k = 0, is obtained by inserting the 
"truncation" © into the exact flow equation 8], for de- 
tails we refer again to 0]. The flow equations for the 
two-body sector, i. e. for the boson propagator param- 
eterized by A^i and m^., are very similar to the SU(3) 
symmetric case (t = \n(k/A)) 
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Since the Yukawa couplings hi are not renormalized. 

d t hi = 0, 
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we can immediately integrate the equations ([3]). The so- 
lution can be found in [3J] . The microscopic values m|j (A) 
(bare couplings) have to be choosen such that the phys- 
ical scattering lengths (at k = 0) between two fermions 
(renormalized couplings) are reproduced correctly. They 
are given by the exchange of the boson field tj>. For ex- 
ample, the scattering length between the fermions 1 and 
2 obeys 
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where all "flowing parameters" are evaluated at the 
macroscopic scale k = and for [i = 0. We use this de- 
scription for the scattering between fermions ip in terms 
of a composite boson field (p also away from the res- 
onance. We emphasize that the field <j> is not related 
to the closed channel Feshbach molecules of the nearby 
resonance. It rather describes an additional "effective 
boson" which may be seen as an auxiliary or Hubbard- 
Stratonovich field, allowing for a simple but effective de- 
scription. For the numeric calculations in this note we 
will use large values of h 2 on the initial scale A. This 
corresponds to pointlike atom-atom interactions in the 
microscopic regime. 

Quite similar to the scattering between fermions ip in 
terms of the bosonic composite state <f> we use a descrip- 
tion of the scattering between fermions ip and bosons <p in 
terms of the trion field \. As an example, a process where 
the fermion tpx and the boson (pi scatter to a fermion ip2 
and a boson <p 2 , is given by a tree level diagram as in Fig. 
[3] For vanishing center-of-mass momentum the effective 
atom-boson coupling reads 
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The flow equations for the three-body sector within our 
approximation are given by the flow of the "mass term" 
for the trion field 
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and the Yukawa-like coupling gi with flow equation 
g 2 h 2 h 1 k 5 (6fc 2 -5/x+5jr) 
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The flow equations for g 2 and g% can be obtained from Eq. 
(|8"|) by permuting the indices 1, 2, 3. For simplicity, we 
neglected in the flow equations and {SJ a contribution 
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that arises from box-diagrams contributing to the atom- 
boson interaction. As described in 0] this term can be 
incorporated into our formalism using scale-dependent 
fields. Also terms of the form tp^ip^^j with i ^ j, that 
are in principle allowed by the symmetries are neglected 
by our approximation in Eq. @. We expect that their 
quantitative influence is sub dominant as it is the case 
for the SU(3) symmetric case (Io| . 

We apply our formalism to 6 Li by choosing the initial 
values of mL at the scale A such that the experimentally 
measured scattering lengths (see Fig. [3]) are reproduced. 
For A t j 1 i(A) = 1, the value of hi parameterizes the mo- 
mentum dependence of the interaction between atoms on 
the microscopic scale. Close to the Feshbach resonance it 
is also connected to the width of the resonance hf ~ A£>. 
We choose here equal and large values for all three species 
hi = /12 = hj, = h. This correspond to pointlike interac- 
tions at the microscopic scale A. Since the precise value 
of h is not known, we use the dependence of our results on 
h as an estimate of their uncertainty. The initial values 
of the couplings m 2 , and gi are parameters in addition to 
the scattering lengths which have to be fixed from experi- 
mental observation. For equal interaction between atoms 
ip and bosons <fi in the UV, the parameter to be fixed is 



A (3) = _ 



g 2 (A) 

m2(A) 
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with g = gi = 92 = .93- Pointlike interactions at the 
microscopic scale may be realized by m^, (A) — > oo. 



state of three atoms X='0i^2V'3- The binding energy per 
atom E of this bound state is given by the chemical po- 
tential \n\ with [i fixed such that m 2 , = 00. To compare 
with the recently performed experimental investigations 
of 6 Li we adapt the initial value A^ 3 -* such that the 

appearance of this bound state corresponds to a magnetic 
field B = 125G, the point where strong three-body losses 
have been observed. Using the same initial value of A^ 3 - 1 
also for other values of the magnetic field, all microcsopic 
parameters are now fixed. We can now proceed to the 
predictions of our model. 

First we find that the bound state of three atoms exists 
in the magnetic field region from B = 125G to B = 
498G. The binding energy per atom E is plotted in the 
lower panel of Fig. [3] We choose here h 2 = 100 , 
as appropriate for 6 Li in the (l,2)-channel close to the 
resonance, while the shaded region corresponds to h 2 6 
(20 a^ 1 , 300 ao 1 ). 

As a second prediction, we present an estimate of the 
three-body loss coefficient A3 that has been measured 
in the experiments by Jochim et al. [4| and O'Hara et 
al. Q. For this purpose it is important to note that 
the fermionic bound state particle x rnight decay into 
states with lower energies. These may be some deeply 
bound molecules not included in our calculation here. 
We first assume that such a loss process does not depend 
strongly on the magnetic field B and therefore work with 
a constant decay width T x for the bound state x- The 
decay width T x appears as an imaginary part of the trion 
propagator when continued to real time 
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FIG. 3: (Color online) Upper panel: Scattering length a\i 
(solid), a,23 (dashed) and 031 (dotted) as a function of the 
magnetic field B for 6 Li. These curves were calculated by 
P. S. Julienne [U and taken from Ref. 0]. 
Lower panel: Binding energy per atom E of the three-body 
bound state x^V'iV^Tfe- The solid line corresponds to the 
initial value h 2 = lOOa^" 1 , while the shaded region gives the 
result in the range h 2 — 20a^ 1 (upper border) to h 2 — SOOa^ 1 
(lower border). 

We solve the flow equations |3j), {4|, ([7|) and ([8|) numer- 
ically. For some range of A*- 3 ) and /j, < we find m 2 , — 
at k = for large enough values of the scattering lengths 
o-i2 1 «23 and 031. This indicates the presence of a bound 
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Instead of working with negative \x chosen such that 
777,2 _ ag ^Qjjg f or t ne computation of the binding en- 
ergy, we now perform an energy shift such that the zero 
energy level corresponds to the open channel and there- 
fore /x = 0. In the region from B = 125G to B = 498G 
the energy gap of the trion is then negative m 2 , < 0. 

The three-body loss coefficient A3 for arbitrary T x is 
obtained as follows. The amplitude to form a trion out 
of three fermions with vanishing momentum and energy 
is given by Xa=i higi/m^. The amplitude for the transi- 
tion from an initial state of three atoms to a final state of 
the trion decay products (cf. Fig [T]) further involves the 
trion propagator that we evaluate in the limit of small 
momentum p 2 — (^iPi) 2 —> 0, and small on-shell atom 
energies uji — p 2 , ui — • u/j — > 0. A thermal distribu- 
tion of the initial momenta will induce some corrections. 
Finally, the loss coefficient involves the unknown vertices 
and phase space factors of the trion decay - for this rea- 
son our computation contains an unknown multiplicative 
factor ck- In terms of p given by Eq. ([1]) we obtain the 
three-body loss coefficient 
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Our result as well as the experimental data points [4j 
are shown in Fig. [4] The agreement between the form of 
the two curves is already quite remarkable. 
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FIG. 4: (Color online) Loss coefficient K3 in dependence on 
the magnetic field B as measured in 0| (dots). The solid 
line is a two-parameter fit of our model to the experimental 
curve. We use here a decay width T x that is independent of 
the magnetic field B. 

We have used three parameters, the location of the 
resonance at Bq = 125G, the overall amplitude ck and 
the decay width T x . They are essentially fixed by the 
peak at Bq = I25G. The extension of the loss rate away 
from the peak involves then no further parameter. 

Our simple prediction involves a rather narrow second 
peak around B\ ss 500G, where the trion energy be- 
comes again degenerate with the open channel, cf. Fig. 
|3l The width of this peak is fixed so far by the assumption 
that the decay width T x is independent of the magnetic 
field. This may be questionable in view of the close- 
by Feshbach resonance and the fact that the trion may 
actually decay into the associated molecule-like bound 
states which have lower energy. We have tested sev- 
eral reasonable approximations, which indeed lead to a 
broadening or even disappearance of the second peak, 
without much effect on the intermediate range of fields 
150G < B < 400G. 

In conclusion, a rather simple trion exchange picture 
describes rather well the observed enhancement of the 
three-body loss coefficient in a range of magnetic fields 
between 100G and 520G. A similar trion dominated 
three-body loss is possible for large B (B > 850G), where 



also a trion bound state with energy below the open chan- 
nel exists. However, the dimer bound states are now 
above the open channel level, such that the trion decay 
may be strongly altered. The role of trion bound states in 
the resonance region is an interesting subject by its own, 
that can be explored by our functional renormalization 
group methods with an extended truncation. 
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